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The basic character of diffusive transport in a magnetised plasma de- 
pends on what kind of transport is modeUed. ExB turbulence under drift 
ordering has special characteristics: it is nearly incompressible, and it can- 
not lead to magnetic flux diffusion if it is electrostatic. The ExB velocity 
is also related to the Poynting energy flux. Under quasineutral dynamics, 
electric fields are not caused by transport of electric charge but by the re- 
quirement that the total current is divergence free. Consequences for well 
constructed computational transport models are discussed in the context of 
a general mean field analysis, which also yields several anomalous transfer 
mechanisms not normally considered by current models. 
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I. Introduction — Transport as Turbulent Diffusion 



Transport of thermal energy and particles in magnetically confined, laboratory plas- 
mas is well known to be anomalous: much larger than transport by coUisional diffusion 
and appearing with different scaling characteristics [1,2]. Modelling of this transport is 
usually done in terms of the two fluid Braginskii equations [3] , but by substituting either 
empirical or theoretical models for the diffusivity coefficients [4]. The transport process 
is assumed by the models to have the same character as the coUisional diffusion process, 
including the prospect that Onsager symmetry should hold [5]. It is important to note 
that the transport process is assumed by such a model to have the same character as the 
coUisional diffusion process. Many of the results of kinetic transport theory rest on the 
assumption that fluctuations in the thermodynamic state variables are small enough to be 
neglected in the energy budget and that they are randomly correlated with each other. 

Small scale, low frequency drift turbulence involving eddies, waves, or vortices of the 
ExB velocity, 

= ;^BxV</., (1) 

where at drift scales the perpendicular electric field is = — V(^, is a fundamentally 
different process from this. Rather than involving dissipation directly, the ExB eddies 
advect the background thermal gradient to produce disturbances in all the thermodynamic 
state variables. The disturbances are then carried with the fiows for as long as the eddies 
last. Then, they are picked up by new eddies and carried further. This does have the nature 
of a diffusive process if the size of the eddies and their lifetime are small and short compared 
to the scales of the background, and the net time-averaged transport does indeed scale with 
the background gradients if the turbulence is local [6]. Moreover, although the dynamics of 
the turbulence is robustly nonlinear, the relative amplitudes of the disturbances are indeed 
small because the plasma is magnetised; the relative amplitudes e^/Te or Pe/Pe-, for the 
electrostatic potential and electron pressure for example, need only be so large as the ratio 
of the local drift scale, ps, to the background gradient scale length, Lp, both defined by 



2 C^M.Te 1 
Ps^ g2g2 Lp=|VlogPe| (2) 



in order to be so nonlinear. Typical measured values of these fiuctuation amplitudes are 
of order one to ten percent, rising towards the edge of a confined plasma due to the 
fact that Lp becomes much smaller than the nominal parallel length scale qR (the field 
line pitch parameter times the major radius of a toroidal plasma) and the perpendicular 
drift dynamics is more able to compete with the dissipation incurred through the parallel 
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electron dynamics [7]. Nevertheless, the strong nature of the coupling mechanisms between 
the state variables (especially Pe and 0) means that the basic properties of this turbulence 
are those of nonlinear drift waves [8,9], under drift ordering — the scale of the turbulence 
is short compared to that of the background gradients and the relative amplitude of the 
fluctuations in the thermodynamics quantities is small [10]. Drift wave turbulence involves 
dynamics on scales which are proportional but not necessarily equal to ps for space and 
Cg/Lp for time, for plasmas which have both the local drift parameter, 5, and the electron 
dynamical beta, /5e, as small parameters. These two parameters and the sound speed, c^, 
are given by 



^=z; = ^ = ^ '' = Mi 

where va is the Alfven velocity. The presence of Tg and reflects coupling between 
electron thermal dynamics and ion inertia. Even with flnite ion temperature Ti ^ T^, the 
coupling between the electron pressure and the electric field keeps ps and Cs/Lp as the 
typical scaling parameters of the disturbances. This coupling and the fact that one of the 
dynamical state variables (0) serves simultaneously as the potential for the electric field 
and the stream function for the (essentially incompressible) flow eddies are the two most 
important properties of drift wave turbulence. These properties and the result that the 
disturbances are not statistically independent of each other have consequences concerning 
the nature of the transport process. Our interest in this article is to explore how these 
properties differ from those of the kinetic diffusion process behind classical (and neoclas- 
sical) transport processes, and what implications this has for well constructed transport 
models. 

II. The Basic Character of Small Scale ExB Flow Transport 

There are a number of fundamental properties of ExB drift turbulence which suggest 
that the transport they cause should be qualitatively different from random thermal trans- 
port. First, if conventional drift ordering holds the velocity is nearly incompressible. If 
the magnetic field is straight then the divergence V ■ v^; vanishes entirely. In a toroidal 
device, the divergence carries the scale of the major radius: 

V • V£; = V ■ -^BxVcj) « -ve • VlogS^ ~ R-^ve (4) 
B 

Since in the confinement zone the pressure scale length, Lp is much smaller, we have 

vs-Vp>pV-Vs if Lp^R (5) 
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Here, we assume that the scale of a significant divergence of the transport flux is Lp rather 
than the much smaUer scale of ExB disturbances, because a locally divergence free state 
for the fluxes should be established on the relatively fast time scale of the turbulence, 
not the transport. By contrast, a diffusive flux of the form tt-cVd = —DVue is always 
compressible, as it can be derived from a potential. 

The second important characteristic is that the ExB velocity is also involved in the 
Poynting energy flux, 

c ^ „ 1 ExB 

with p the mass density. Due to the close force balance for small scale disturbances in 
drift dynamics, there is a cancellation between the Poynting flux and one factor of pvg 
in the total thermal energy flux. This has the consequence that alone among transport 
mechanisms, for ExB drift turbulence the advective part of the thermal energy flux appears 
with a factor of 3/2, not 5/2, times the temperature times the particle flux [11]. More 
speciflcally, its energetic coupling to the thermal reservoir is mediated through its small, 
quasistatic divergence, so that the pressure does little work on elemental volumes of the 
fluid. 

A third characteristic of ExB turbulent transport also arises from the connection be- 
tween the flow stream function and the electric fleld: anomalous energy exchange channels 
between the electron and ion thermal energies. The basic mechanism underlying drift 
waves is the coupling between electron pressure and electrostatic disturbances through the 
parallel current: Pe ^ J\\ ^ 4'- The second step in this is the Alfvenic coupling Jy ^ (p 
which is present in MHD. The flrst step is the adiabatic coupling pe ^ Jy, part of which 
is the adiabatic response in the Ohm's law, by which pressure disturbances launch shear 
Alfven waves along the background magnetic field due to the force balance for electrons 
parallel to the fleld. This adiabatic coupling is responsible for a high degree of correlation 
between the electrostatic potential and electron pressure disturbances, and it provides the 
means by which the free energy in the background gradient is passed to the ExB turbu- 
lence [6,9]. It can only be addressed by a proper two fluid model which treats electron 
and ion dynamics separately but self consistently. Since the ExB drift kinetic energy is 
essentially in the ions, and the conservation of charge keeps the dynamics quasineutral, a 
parallel current divergence in the Alfven dynamics is balanced by a polarisation current 
divergence in the ions, and this can do work on the ion fluid. Through the disturbances, 
electron thermal energy taken out of the background gradient can be deposited in the ions. 
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A fourth characteristic of ExB flows is that they cannot cause magnetic flux diffusion 
unless there is a significant amount of reconnection of magnetic field lines in the presence 
of a current gradient (that is, the drive mechanism should be this current gradient, rather 
than the thermal gradients as for drift turbulence). The magnetic flux occurs through 
the inductive part of the electric field in the force balance for electrons, that is, dA\\/dt. 
Neglecting dissipation and electron inertia the Ohm's law is 

In an average over a closed flux surface the right hand side vanishes unless there are 
appreciable magnetic disturbances. When the disturbances follow from drift dynamics the 
relative phase shifts are such that magnetic transport effects tend to cancel out of the 
flux surface averages. Only when the current gradient is available as an energy source do 
these processes lead to appreciable transport, even if the thermal transport is robustly 
anomalous. By contrast, random thermal diffusion transports magnetic flux by the same 
mechanism as for particles and energy, and the flux diffusivity, ?7||C^/47r, is even larger than 
the fluid diffusivities, proportional to D^, = pl^ei by a factor of f3~^. 

Two further considerations involve the development of large scale ExB flows within the 
flux surfaces (we assume that the equilibrium is quiescent enough that there is no bulk flow 
across flux surfaces beyond the existence of the turbulence). For dynamics at drift scales, 
the disturbances as well as the background are deeply quasineutral. Charge differences are 
neutralised on time and space scales which are effectively arbitrarily small, so that while 
there can be nonvanishing electric fields, there is no significant charge density. In this 
context it is inappropriate to speak of a radial electric field generated by the divergence 
of a radial current. Rather, we have the constraint that the total current is divergence 
free, so that various pieces of the current which might have a divergence are balanced by a 
divergence in the polarisation current of the ions. This is relatively well known in a general 
sense, but we will illustrate it here specifically for the problem of mean fiow generation. If 
source or loss mechanisms for either of the charge species are considered, these should be 
thought of as a torque on the ExB vorticity rather than a charge generation[12]. 

Finally, anomalous transport in the ExB vorticity results fundamentally from anoma- 
lous viscosity, either the fluid Reynolds stress or some correction arising from gyroviscosity. 
Under flux surface and short time scale averaging, it is a matter of the total radial polari- 
sation drift (linear plus nonlinear) necessarily vanishing. An anomalous resistivity cannot 
lead to anomalous transport of ExB flows or, equivalently, generation of a radial electric 
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field. This is a direct consequence of the fact that resistivity is a momentum conserving 
friction between electrons and ions. 

In the rest of this article, we address the general equations of transport by small scale 
drift turbulence via mean field theory, and then underscore the above points via analysis 
of these equations. 

III. Global Fluid Drift Equations and Energy Conservation 

The general derivation of fiuid drift equations for large scale motion proceeds similarly 
to the local treatments more familiar from the turbulence studies [8,9]. The the difference 
is mainly that we must keep track of the variability of the coefficients which depend on 
densities or temperatures. The density is the most important of these. It carries the 
consequence that the ion polarisation drift must be retained in advection if exact energy 
conservation is to hold. This result follows from the fact that proper consideration of the 
electron force balance (the adiabatic response) disallows the use of the more familiar MHD 
ordering in the derivation of global equations [13,14,15]. What we must do is to generalise 
sufficiently to retain a consistent superset of both approaches. 

We begin with a neutral, single-component plasma with ions of mass Mi and charge 
state Z and electrons of mass mg. For this treatment we will neglect electron inertia but 
keep me in combination with z/g in the formulae for the dissipative fluxes. The densities 
satisfy = Zrii. We assume that the local drift parameter, 5, and the electron dynamical 
beta, /3e, defined in the Introduction, are small. In contrast to the familiar situation with 
drift equation models, no assumption concerning fiute mode ordering in the derivatives 
is made during the derivation of the equations themselves and their energy theorem con- 
cerning spatial scales, even for B. The more familiar assumption of /cy ^ or for a 
toroidally confined plasma Lp <^ R, is reserved for evaluating the properties of the trans- 
port considering that the disturbances involved in the turbulence do indeed satisfy local 
drift ordering. 

The drift approximation basically states that perpendicular compressional wave dy- 
namics is arbitrarily fast compared to the physics of interest and is therefore neglected 
in favor of a quasistatically evolving force balance perpendicular to the magnetic field 
lines. We therefore have a quasistatic perpendicular electric field, with magnetic induction 
appearing only in the parallel component, 

E^ = -Vi0 £;|| = ___J_V||<^ (8) 
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Correspondingly, the total magnetic field, B^, is given by the background, plus the principal 
disturbance arising from ^y, plus a small compressional piece which evolves quasistatically, 
plus an even smaller correction needed to ensure that V • = 0, 



B ~ B 

B^ = B - - X V^ii +B-- V±XB (9) 

respectively. The symbols B and B give the equilibrium magnetic field and its magnitude, 
while b is reserved for the combination of B and the piece involving A\\, given by 

B^ = -|xV^II (10) 
so that in describing parallel components we use 

b=l(B + B^) V||=b-V (11) 

Perpendicular components are described by the background field only, 

BxV o , , 

Vi = -Bx-^ Vi = V-Vi (12) 

so we must note that with these definitions we cannot write V = V_l + bVy, for example. 
The Ampere's law involving Ay is written with V^, 

Att 

Vi^ll = --Jll (13) 

The compressional part of the magnetic field, involving B, is written for generality but 
not used in a /9e <^ 1 treatment (it can be used to evaluate the Poynting energy flux in 
the energy theorem by assuming that of the sum B^ + 8np is zero). The divergence 
correction, involving xb^ is never used; it is written only to preserve V ■ B^ = 0, and in 
practice it is always negligibly small. 

The quasistatic force balance allows us to write the drift velocities and heat fluxes 
perpendicular to B in terms of the state variables. For electrons the diamagnetic fluxes 
are 



and for ions 



The ExB drift for both species and the diamagnetic current are given by 

V£;=-^BxV0 J* = nee(u* - V*) = -^BxVp (16) 

with p = Pe + Pi- We form the fowest order perpendicular fluid velocities resulting from 
these drifts, 

U_L=Vs + U* V_L=Vs + V* (17) 

In the parallel dynamics the dependent variables are u\\ and Jy, with v\\ defined as 

"II = n-^, (18) 

The parallel current is actually an auxiliary variable obtained from Ay through the drift 
form of Ampere's law, Eq. (13) above, with replaced by V^. 

For the ions we keep the polarisation drift, Up, which is required in the current balances 
and also in the contributions of the ion pressure to energy conservation. The polarisation 
drift is obtained by inserting the lowest order ion velocity, Uj^ + Mub, into the ion inertia 
terms (including all components of the stress tensor), and solving anew for the total ion 
velocity, 

u = u_L + tt||b + Up (19) 
For the electrons we neglect polarisation along with mg, so that their total velocity is 

v = v_L+V||b (20) 

noting that f y is given in Eq. (18). The polarisation drift is given by 

+ -|^Bx(V-n.) (21) 

where we keep the gyroviscosity explicitly instead of just using the part which cancels 
■ V, in order to preserve its symmetry. We can express the polarisation drift in terms 
of a potential, 

following from the fact that it is primarily proportional to gradients of scalars. In contrast 
to the local treatments, the polarisation drift does indeed appear explicitly in the inertia, 
but only in the advective derivative. On the other hand, the inertial velocity itself is the 
lowest order version, whose perpendicular component is Uj^. Consequently, we construct 
the components of 11* using only the lowest order velocities, u_l + tt||b, as we will need to 



c f d 
UeSUp = UiMi— Bx I — + u • V 



V0 + 



Ue.e 
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use the fact that 11* : Vu vanishes, for whichever vector u is used to construct 11* (if this is 
done in detail it must be extended to include qj/\ in 11* to get the gyroviscous cancellation 
[13] to work properly when the temperature and density gradients are equally steep [16]). 
We note here that everything having to do with ion inertia in the perpendicular force 
balance is encompassed within the polarisation drift, including not only the advection of 
vorticity but also Reynolds stress and gyroviscosity. 

The equations under the drift approximation are formed by substituting the above 
expressions for the velocities into the rest of the fluid equations. The evolution of ^ is 
controlled by the polarisation equation, which is formed by subtracting the densities so as 
to form a continuity equation for electric charge. The charge density itself is set to zero, 
leaving V ■ J = 0. Expressing the total current as 



J = Jllb + J* + Jr 



(23) 



with the polarisation current given by 3p = UeCUp, we find 

-V • Ueeup = V • (J||b) + V • BxV (pe+Pi) 
Where we need to evaluate Up in the ion equations, we may use. 



S2 



(24) 



(25) 



which we solve for x- The evolution of and Jy is controlled by the Ohm's law, which 
is just the parallel electron force balance including dissipation. The rest of the dependent 
variables in the set {(f),ne, A\\,pe,u\\,pi} have their own fluid equations, in which we use 
the same procedure as in the polarisation equation. The resulting model is given by 



+ — Bx(v-n*)^ = 



d_3 
di2 



due 
dt 



+ V 



V • (J||b) + V • Bx V (pe +Pi) (26) 



Jii 



1 dA\\ _ V\\Pe 



B 



b--^BxV^ 
ripe J B'^ e 



c dt 



UeC 



-Vn6-Rp 



(27) 
(28) 



Pe + V 



lPe-E + lpe[U\l-^^^ 



b + 9e||b--^BxV 



5 peTe 

2 e 
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= W||V||Pe 



V^Pe-Pe"^ ■VE + J\\Re^ (29) 



„.M. (I ^ u . v) „„ + b . (V . n.) ^ -V„ (p. + V . (.„bV„„„) (30) 



_a 3 
dt2 



5 
2^ 



PiV^; + -PiU\\h + gi||b + -PiUp + — Bx V 



5 piTi 
2 e 



_ __ I __ I Z 

W||V||]3i + Up- Vpi-]3iV-V£;+;U|| |V||W||| (31) 



where except for V • the magnetic divergences are written exphcitly. The dissipative 
momentum transfer, 
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a I ^ e 
J\\-\ aJ\\ + —qe 



and parallel heat fluxes, 



Qell +a J\\ = -Kene — V\\T, 



Vp 



qui = -KiUi—VnTi 



(32) 



(33) 



are given by the Braginskii model, with the resistivity and electron and ion thermal veloc- 
ities given by 



me^e 
= ^ 2 



rUe 



1/2 



K- = ( — 

Mi 



1/2 



(34) 



respectively, with numerical constants for Z = 1 



77 = 0.51 



a 



0.71 



3.2 



3.9 



(35) 



and the collision frequencies given by the inverses of the Braginskii collision times [3] . 

We note that this is an approximate treatment, but as we will now show, it conserves 
global energy. We first make a few remarks. First, the reason for including Up in the 
advective derivative u • V is conservation of the inertial energy incorporated in u_l. We 
must manipulate the polarisation equation such that the density is brought under both the 
time derivative and under a resulting total divergence. Since it starts out as a multiplier 
on d/dt as well as on the advecting velocity, the continuity demands require that we keep 
Up in riiU • V if we keep it in V • njU. This chain of demands starts at the presence of VyPe 
in the Ohm's law, continues with Vy Jy in the electron pressure equation and hence the 
density equation, and finally with the fact that the divergences of Jy and Up go together 
in the continuity. The ordering and conservation are preserved by keeping Up in u • V, but 
not in the actual inertial velocity, u_l. Second, although the perpendicular and parallel 
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components of vectors are set up differently, they are kept separate in all computations 
so that there is no violation of energy conservation. Finally, one sees that part of the 
drift ordering is expressed in the fact that the equations are written only for the lowest 
order components of the ion velocity and magnetic field, but the operations involving those 
velocities and gradients include all components in order to conserve energy. 

The energy theorem is formed by multiplying Eqs. (26,30) for the vorticity and parallel 
velocity by —cj) and u\\, and Eq. (28) for the magnetic potential by Jy, respectively. The 
pressures in Eqs. (29,31) represent the thermal energy. The total energy density is 



= j (fx ^niM^i^ 



U= / d^X I UiMi^^ + mMi^ + -pe + + ) (36) 



We will call these pieces the drift energy, sound wave energy, electron and ion thermal 
energy, and magnetic energy, respectively. 

The first consideration is how polarisation works under the drift approximation. We 
multiply the polarisation equation by (j), to obtain 

0V • J = V • 0J - J • V0 = (37) 

where J includes all three pieces, Jyb + J* + Jp. The form of this that we actually use is 

-Jp- V0 + V-(/)J- J* • V0 = J||V||(/) (38) 

where the first piece will give part of the drift energy, the next two pieces will be combined 
into a total flux plus a transfer term, and the piece on the right side is already a transfer 
term. We manipulate 

J* • V0 = (J^ Bx Vp) • = - (-^ Bx V0) • Vp = -ve ■ Vp (39) 

where p = Pe + Pi is the total pressure, and then combine with the total divergence to 
write 

-Jp- V(/) + V- ((/)J+pve) = J||V||(/)+pV ■ (40) 
We note that in general the divergence of is equivalent to part of the Poynting fliix, 

V >J = V • (^0VxBt) = -V • (^V(/>xBi) (41) 

dropping a divergence of a curl. The near cancellation with p\E is called the Poynting 
cancellation, 

V- (0Ji+pvs) «iO (42) 
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affecting since is also perpendicular to B. This cancellation results in the fact 
that the ExB thermal transport eventually appears with a factor of 3/2, so that the ExB 
divergence, not an ExB advection, is the process by which transfer with the thermal 
reservoir occurs. We will show this in Section V, together with the other characteristics of 
ExB transport. 

The part of the Poynting flux involving Jy moves parallel to B and is the process by 
which the MHD part of shear Alfven waves move field energy along B: 

V • (0J||b) ^ -V ■ (^V0xBi) (43) 

neglecting the action of the curl on b itself. This is generally comparable to or larger than 
V ■ (pet'||b), reflecting the properties of shear Alfven dynamics under general two fluid 
ordering. 

The pressure piece in the drift energy is found by subtracting Up ■ Vpi from both sides 
of Eq. (40), to find 

-Jp • {V(f> + + V • (0J + pvs) = ^11 V||0 + {pe + Pi) V • V£ - Up • Vpi (44) 

Finally, substituting in for Jp and noting how the combination in parentheses gives rise to 
u_L, we have 



dt\ 2 ' 



UiMi ^ U + U_L • + (j)3 -\-pVE 



= : Vui + J\\V\\(t)+ ipe + Pi) V • - Up • Vpi (45) 

in which the transport terms are the ones under the divergence operator on the left side, 
and the transfer terms are those on the right side. We note that the reason that u must 
include all components of the ion velocity is to bring the factor of under both the partial 
time derivative and the divergence operator simultaneously. In other words, whatever we 
keep in the velocity divergence in the density equation must also be kept in the advection 
in the polarisation equation, and in all the ion fluid equations as well. 

We now find where the energy transferred out of this drift energy goes. The next 
subtlety is the magnetic energy. We multiply Eq. (28) by J\\ and manipulate the divergence 
to obtain 



~ |V.^„r + V . (Vx^i,^^) = iv,p, - J„V„^ - (46) 
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The terms on the right side are the adiabatic and Alfvenic transfer effects and resistive 
dissipation, which function the same way as in drift wave turbulence [6]. The divergence 
term is the inductive part of the Poynting energy flux, as we can see by evaluating the 
projection operations implicit in the V_l operator. 



V_L^|| TT^ = V 

'Uttc dt / 47r 



— xVAi 



(47) 



cB dt \B 

Substituting B_l for A\\, this gives the conservation law for the perturbed magnetic energy 
d IB^I^ _ c /I, dAu ^ \ J|| 



dt 8tt 



where we note that the inductive electric field appears with b since the term it is crossed 
into is Bj^. The inductive part of the Poynting flux involving Bj^ is in this equation, 
transporting magnetic energy, while the static part is in Eq. (43), transporting ExB en- 
ergy. Both are small, however, except for the part propagating along B, which carries the 
Alfven wave energy but whose contribution to transport across magnetic flux surfaces is 
nevertheless negligible. 

The sound wave energy is found by multiplying Eq. (30) by wy, noting that here also 
the need to put the factor of Ui under the time derivative and the divergence requires 
keeping Up in the advection terms. We obtain 

d f f ul \ 

— ImMi^ I + V - I niMi^u + W||b-n* - it||b//||V||W|| I = 

= n*:V (M||b) -wiiVii (pe+Pi) - l^w |V||M||| (49) 

noting that the gyroviscosity term cancels properly with the one in Eq. (45) because 
n* : V(u^ + W||b) vanishes. Gyroviscosity thereby represents a diamagnetic momentum 
flux which acts to transfer energy between parallel and perpendicular fluid motion. 

The thermal energy comes next. These are simply the pressure equations, with some of 
the advection terms evaluated piece by piece. For the electrons we have Eq. (29), in which 
the transfer effects with other equations are now obvious: adiabatic transfer (Jy VyPe) and 
resistive frictional heating (J||i?ei) with the magnetic energy, acoustic coupling (tt||V||Pe) 
with the sound waves, and the ExB divergence (peV • v^;) with the drift energy through 
polarisation. For the ions we have Eq. (31), in which the transfer effects are thermal 
coupling (■U||V||Pi) and viscous heating (the term) with the sound waves, and the ExB 
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divergence (p^V • v^;) through polarisation and the polarisation advection (up ■ Vp^) with 
the drift energy. 

The total energy theorem is then given by 



3 

^{Pe+Pi)vE + (ge|| b + 



+ o + Pev\\ + PiU\\) b + -piUp + BxV-^--2BxV-^ + 



|u±|^ + «M ^ ^ ^ . c 1 aAii 

+ niMi ^u + n, • (ui + W||b) + 0J±+pv£ - b— /xB^ 

2 V n / 4:77 C Ot 



(50) 



The transport effects are written in order of their usual importance; in practical situ- 
ations only the terms in the first line need be kept. In situations where transport by 
Pfirsch-Schluter currents (generally caused by nonzero magnetic divergences involving the 
equilibrium pressure) is relevant, the second line would have to be kept as a unit, since the 
divergences of all the currents go together. The diamagnetic fluxes are combined into the 
terms involving the drift operator. The terms in the last line are negligible unless there 
are transonic or supersonic flows. 

The dominant transport effect across the magnetic flux surfaces is the ExB advection, 
(3/2)peV£;, which appears with the factor of 3/2 due to the Poynting cancellation, Eq. (42). 
The diamagnetic fluxes drop to the second line due to the well known diamagnetic can- 
cellation in the pressure equations [17]. For the parallel transport, the principal effects 
are the electron and ion heat fluxes, B^^{qe\\ + qi\\)'B±. We also have the ion advection, 
{5/2B)piU^\'B±, appearing with the factor of 5/2 due to the compressibility of the parallel 
flows. For the electrons the part of this effect involving is small because Ay is the stream 
function for B^ and it is also related to Jy through Ampere's law, but here we also have 
the acoustic advective effect, {5/2B)peU\\B±. These parallel advective effects are small, 
however, since the sound wave transit time is generally much longer than the time scale of 
the turbulence. 

To summarise, we have a set of drift equations following the drift ordering as to 
magnetic compression but allowing for an arbitrary scale of motion and with it, interactions 
with the thermal gradient of the background. The conservation of total energy is exact, 
but at the price that the total velocity must be kept in advection. The extra piece, that is, 
the polarisation drift, is given in terms of the dependent variables through a constitutive 
relation, Eq. (25). 
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IV. A One Dimensional Mean Field Model of Transport by ExB Turbulence 



In this section we construct a mean field transport model including both perpendicular 
and parallel flows. All dependent variables are split into profile quantities which are flux 
functions (dependence only on the flux surface label coordinate), and disturbances which 
have arbitrary coordinate dependence but are zero in ensemble average. The proflles are 
assumed to be slowly varying in time and space compared to the disturbances, so the 
ensemble average can be understood as one over an interval of time short compared to the 
transport time of a given radial region but long compared to the fluctuation time scale 
of the turbulence. The energy content of the disturbances is assumed to be negligible, 
so that their energy equations become statements that their net transfer effects are all 
in quasistatic balance. Transport of wave energy across flux surfaces (disturbance energy 
content times v^;) is kept until the disturbance energy equations are evaluated, in order 
to show that it is small. The disturbances form nonvanishing quadratic transport and 
transfer quantities in the mean field equations; these quantities control the evolution of 
the profiles. 

We retain both poloidal and toroidal rotation, fully described by the two state vari- 
ables, (j) and Pi, and the parallel fiow u\\. In contrast with more usual treatments, we keep 
the partial time derivatives in the polarisation equation, and hence also the energy content 
of the drift flows. Both flows are affected by Reynolds stress and by neoclassical friction, 
the latter acting only on the poloidal component of the total flow (recall that parallel is not 
purely toroidal and perpendicular is not purely poloidal). The parallel flow disturbances 
are also dissipated by parallel viscosity. We also retain magnetic flux diffusion, keeping 
resistive dissipation for both profiles and disturbances. 

The equations for the profiles are those derived in Section III, simplified by the symme- 
try assumptions, but with the addition of the ensemble averaged transport effects. For the 
disturbances we mainly have to declare what is kept in order to satisfy energy conservation, 
since the role of the disturbances is to provide transfer channels. 

To proceed with the transfer evaluations we need the density equation. This model 
follows transport of electrons, but the ion density is needed for the factors of UiMi in front 
of the partial time derivatives in the fiow equations to conserve energy properly. For the 
electrons we have 

^ + V • (nev^) = (51) 
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neglecting compression effects. With this simphfication, we have equivalently the same 
equation for the ions 

^ + V • {uiVE) = (52) 

We write both these equations to emphasise the fact that neglect of the polarisation com- 
pression and the neglect of particle transport by the parallel current in general go together. 
It is important here to note the assumption that V ■ ?ijV^ can be neglected, since the 
density- weighted average of v^; should vanish in order to have transport not go by direct 
flow but by statistical average. This helps in the manipulation of the Reynolds stresses. 

The parallel flow proflle equation is 

riiMi^ + Mi (uiVE) ■ Vw|| + V • UiMi {u\\Ve) + niMiVdU%^ = (53) 

in which we have used the ability to neglect V • riiVE before taking the ensemble average. 
The last term is the parallel component of the poloidal friction, with = (u± + ^||t)) ■ 
the total poloidal flow component, = h ■ V9 the poloidal component of the magnetic 
unit vector, z/^ the neoclassical friction coefflcient, and given by the profiles of (f) and 
Pi as in Section III. The third term is the "conductive" parallel momentum transport, as 
we can see by using Eq. (52) to form the parallel momentum equation, 

— [niMiU\\) + V ■ (ivh {u.we) W|| + n,M, (w||Ve> ) + niMiUdU%^ = (54) 

The first term under the divergence is transport via the anomalous particle flux, and the 
second is due to correlations between the parallel flow and ExB disturbances, which is also 
the perp-parallel Reynolds stress. The latter term is called the "conductive" part of the 
transport. The transport operating through the particle flux is called the "convective" 
part. We simply generalise this concept, which is standard for thermal energy, to the 
other conserved quantities. The last term transfers momentum with the magnetic fleld; 
the momentum of Uj^ is not followed, consistent with the drift approximation in which the 
background magnetic fleld is an anchor. 

For the parallel flow disturbances we have 

du\\ _ _ _ , _ _ _ o_ , , 

UiMi-^ + UiMiWE ■ Vu|| -I- Mi {niVE) ■ Vm|| -|- UiMiVE ■ Vu|| = -Vyp -|- V||U|| (55) 

where p = Pe + Pi is the total pressure disturbance. The second and third terms are 
nonlinear advection, which is required in light of Eq. (52) by the need to get the factor 
of UiMi under the time derivative. The fourth term is the drive term arising from the 
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background flow gradient. It conserves energy against the perp-paraUel Reynolds stress in 
the profile equation. The energy equations are given by 




+ V • I Mi {mWE) -j- + TliM, Vs) W|| I = ^^^^ 



TiiMi {u\\Ve) • Vtt|| - niMiVdU^ u\\V 



e 



for the profile, and 




= -n^Mj (w||V£;> • Vw|| - (w||V||p) -//|| (|V||S|| 



(57) 

for the disturbances, in both of which Eq. (52) times the relevant factors of squared velocity 
(specific energy) is used. This also uses V-njV^; = to treat the conventional nonlinearity, 
and in this and the other disturbance energy equations the resulting term in the energy 
fiux is a third order nonlinearity which represents the ExB transport of wave energy in 
the turbulence. We carry these terms here for illustration, but since they are down by two 
orders of S they are always small in practice (and are measured to be small in turbulence 
computations). 

The energy transfer via the Reynolds stress is apparent in Eq. (57), as are the parallel 
fiow transfer efi'ects we saw in Section III. These four Eqs. (53,54,56,57) together show why 
the particle fiux has to be kept in advection; otherwise we would have problems with si- 
multaneous conservation of energy and momentum. The same is true for the perpendicular 
fiow, which we treat next. 

The perpendicular fiow is determined by the scalar fields (f) and pi, and since the ion 
pressure has its own equation this one becomes the equation for (p. For the profiles we 
have 



= V ■ Bx (Mi (niVs) • Vui + V • n^M^ (v^vs) + mMiViU^VQ^ 

in which the only effects are nonlinear forcing and frictional damping, as with the tt|| 
profile. As in Eq. (53), we rewrite the • Vv^; term as a Reynolds stress. The more 
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familiar linear terms and the conventional vorticity nonlinearity are all in the equation for 
the disturbances, 



neglecting the frictional damping due to its slow time scale. We are keeping only in 
the vorticity disturbances, as the ExB contributions to the Reynolds stress are dominant. 
The first term on the right side illustrates ExB shearing but will drop out of the energy 
equation due to symmetry. The next term is the conventional nonlinearity and will give 
rise to the ExB wave energy transport. The third term is the advection of the flow profile 
by flow disturbances, and is the term which conserves energy against the Reynolds stress 
in Eq. (58). 

The energy equation for the flow profile is constructed as in Section III: we multiply 
Eq. (58) by —cj) and manipulate the divergences (equivalent to integrations by parts under 
an integral), also subtracting (Jp/nge) • Vpi from both sides, the left side completing the 
specific energy (squared velocity) and the right side providing the polarisation transfer 
with the ion pressure. Eq. (52) is used to bring the factors of rii under the partial time 
derivatives. For the profile we have 



S2 di 




V.-^Bx(n, 



MiVE ■ Vve + UiMiVE ■ Vvs + UiMiVE ■ Vu± ) (59) 





UiMi^j + V • \^(j)Jp + Mi {tuve) -j- + riiMi (veVe) ■ u^j = 

= UiMi (vbVe;) : VU_L - TliMiVdu'^ u'^^ - Jp ■ 




(60) 



UeC 



where tt' 



u_L • V^, and 



+ Bx (^Mi {ni^rE) • + V • n^Mj (v^v^;) + mMiUdU^Ve 
is the polarisation current, and for the disturbances we have 




(61) 





(62) 
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where u]_ = v|. since we are keeping only ExB vorticity in the disturbances. We wee that 
the Reynolds stress transfers within the flow energy between profile and disturbances, and 
the disturbances then transfer to the pressures (pV ■ Ve) and the magnetic disturbances 
(J|| V||0). The wave energy transport is small by 0{5^) compared to either (pv^) or (0J_l). 
The cross flux surface Alfven wave flux ((f)J\\h±) is also small by 0(5^). 



The electron and ion pressure equations are somewhat easier, since they are already 
written in conservative form. We keep the most important quadratic fliixes, plus all the 
transfer effects of the disturbances. For the electrons we have 

+ ^^^^^^ + (^^11^^) + 1^^ ("11^^)) = ^g3^ 

= (^11 V||Pe> - — ( J|| V||Pe) - (PeV • V^) + V\\jf\ + {j\\R^ 



neglecting the magnetic flutter transport involving Jy . For the ions we have 

Sdpi ^ [?, , /- r \ 5 ~ 
2^ + ( 2 ^^^^^ + \?i||b±^ + -^Pi {u\\h^ 



= («||V||Pi> - (piV • Ve) + (up • Vpi) + (64) 



+ mMiVdU^u'' + — • Vpi + (| Vy-u 

correspondingly neglecting the transport due to Up. This is the form of the pressure 
equations including magnetic flutter transport; formally the part of that due to u\\ goes 
together with the transfers adding up to (tt||V||p), but while the transfer effects remain, 
the transport effects are overshadowed by the contributions to (pv^)) so we drop flutter 
transport at this stage. Since Reynolds stress is not a proper viscosity, it does not a priori 
appear in the ion pressure equation. It represents a transfer effect with the disturbances, 
which in turn transfer the energy via other channels with the ion pressure, most notably 
through sound waves and through some of the electron channels, as we will see. 

The magnetic fleld evolves through fliix diffusion, 

.-l=E,-r„J, (65) 

balanced by the loop voltage E^, here also neglecting magnetic flutter transport (always 
small here unless driven by reconnection, i.e., a current gradient). The energy in the 
disturbances evolves according to 



d_l_ 

dtSTT 



B 




J||V||Pe>-(J||V||0)-(J||i?ei) (66) 
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Figure 1. Energy transfer diagram for the disturbances in drift wave turbulence. 
Large and small scales are indentified as low and high k±, respectively. The 
main state variables are (p and pe and pi , the latter two drawn collectively as p. 
The transfer channel ^ J|| ^ Pe maintains drift wave mode structure, while 
toroidal compression p ^ (p represents interchange forcing, especially impor- 
tant for ITG turbulence since the channel through Jy does not directly affect 
the ions. The weaker transfer channel pe u\\ <-> pi representing the sound 
waves is not shown. Transfer between scales proceeds through the polarisation 
nonlinearity for (j) and the ExB advection nonlinearity for p. The principal sink 
is diffusive mixing of p out of the energy producing spectral range by ExB ed- 
dies, and the resistive sink is also shown. The source is the background pressure 
gradient. 



keeping Rei general as it also contains parallel gradients of disturbances (cf. Eq. 32). 

The key question is what to do with the equations for the disturbed energies in a 
transport model. A reasonable treatment is to note that over an ensemble average the 
partial time derivatives are zero. The energy content of the disturbances is to be neglected 
in any case, and since the wave energy transport terms are all small by 0{S^) we drop 
them as well. These equations then merely become consistency relations among all the 
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transfer effects. Eqs. (57,62,66) respectively become 

(w||V||Pi> + //|| (|V||S||n = -niMj(w||VB> • Vw|| - (w||V||Pe> (67) 



(PiV • ve) = -riiMi {vE^rE):Vui_ + ( J||V||0^ - (peV • v^;) (68) 

•^l^^^) " (•^l^ll^^) = - (•^l^ll^) ^^^^ 
which will be used to evaluate their the left sides in the equations for pi, pi, and pe, 
respectively. In contrast to the turbulence, whose nonlinear character is sensitive to the 
average size of the transfer effects {i.e., the standard deviation) [7], the profiles sense only 
the average transfer and transport (i.e., the mean). Some of the transfer effects cancel 
in this fashion, most notably the adiabatic transfer mechanism (^J^^V ^^Pe"^ ■ The transfer 
dynamics in ExB turbulence is depicted in Fig. 1, with the nonlinear drift wave and 
interchange mechanisms highlighted. All of the arrows are bi-directional, due to the fact 
that the mixing and scattering tendencies in the nonlinearities, which transfer between 
scales of motion, impose their quasi-random character on the linear mechanisms which 
transfer between the state variables Pe, Pi, and 0, in some cases with the flux variables Jy 
and as intermediaries. 

The transport model consists of equations for the usual three thermodynamic state 
variables, plus two more for the flow quantities, plus one more for the magnetic fliix. It is 
written as follows: 

^ + v-r = o (70) 

3 dp. 



2 dt 



+ V-Qe = Tei+ 7711 J,f (71) 



3 dv 3 

+ V • Qi = -Te^ - : Vui - R£|| • Vw|| + n,M,i^dU^u^ + — -Vpi (72) 

^n^^^ (^^^^Y^^ ■ ^Bx[MiT -Vu^ + V ■YLE + niM.Vdu've) (73) 
niMi^ + M^r • V«|| + V • R^ii + niMiVdU%^ = (74) 

The parallel current is given by Ampere's law. 



J\\ = -^^U\\ (76) 
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The following transport quantities are defined as ensemble averages over the turbulence: 

3 3 

r = (nevs) Q« = 2 ^^^^^^ ^'^2 ^^'^^^ ^^^^ 

R^; = UiMi (ve^^e) R^ii = niMi (wyvs) (78) 

Te^ = (S||V||Pe> - (PeV ■ V^) - V||^) (79) 

giving the particle flux, electron and ion heat fluxes, the perp-perp and perp-parallel 
Reynolds stresses, and the anomalous heat transfer, respectively. 




A 



loop voltage 

Figure 2. Energy transfer diagram for the profile quantities in a transport 
model. The five quantities representing energy in Eqs. (71-74) are shown, 
with the transfer channels representing the disturbances. The three pathways 
connecting pe to pi represent the three terms in Tgi in Eq. (79). The source is 
the loop voltage. The only sink is transport through the outer boundary. 
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The perpendicular flow and polarisation current are given by 

u. = ^Bx(v^+^) (80) 

Jp = Bx (n,M,^ + MiT ■ Vu^ + V • + UiMiUau^Ve'^ (81) 

and poloidal contravariant components arc denoted by the superscript 9. As explicit 
dissipation mechanisms we have the neoclassical friction Ud and the parallel resistivity 
?7||. The only source in the model as written is the loop voltage E^. 

The most important result of this model is the set of transfer channels between each 
of the pressures and the perpendicular and parallel flows, which are depicted in Fig. 2. 
This includes the appearance of anomalous transfer effects between the electron and ion 
temperatures, given by Tgi as defined in Eq. (79) . This anomalous transfer is an essential 
result of drift wave physics in the turbulence. The electron thermal gradient drives ExB 
motion through the adiabatic and Alfvenic couplings, Pe ^ Jy (f). Dissipation takes 
place in both the electron (</||) and ion {u\\) channels, which is here represented in terms 
of the coupling between the pressures. More generally, whichever species has the larger 
temperature contributes more to the drive of the turbulence, relaxing the corresponding 
gradient the more efficiently. Rather than a monotonic exchange effect as would be man- 
dated by dissipative kinetic theory, the anomalous transfer is more a tendency towards 
equipartition given a statistical system. All three of the terms which comprise T^i in 
Eq. (79) are involved, but (^Jy V||(/>^ is the strongest for drift wave turbulence, {pe'V ■ ve) 
is strongest for interchange turbulence, and VyPe) is strongest for ITG turbulence (the 
pure toroidal ITG transfer goes through (p^V ■ v^;), i.e., directly from ion flows to the 
ion pressure, not between species). For strictly cold ions in the absence of flows, the bal- 
ances among the disturbances (Eqs. 67-69) would indicate that Jy V||(/)^ would balance 
(PeV • v^;) and that (u\^V\\Pe) would vanish and hence that T^i would be zero on average. 
However, the presence of both flows as reserviors allows a flnite T^i to exist even for = 0. 
Although this probably depends in detail on the existence of direct dissipation in the ions 
(i.e., viscosity), and therefore a finite T^, it is nevertheless important to the physics that 
these transfer channels exist and that they can be significantly stronger that the classical, 
coUisional transfer usually considered. It is important to note that the drift wave anoma- 
lous heat transfer, (^J\^V \\(/)^ , representing compression of the polarisation drift in the ions 
and parallel currents in the electrons, can be significant even for Ti — 0. 

The other significant transfer effects are those involving the Reynolds stresses. The 
one most considered is R^;, acting on the perpendicular flow. But in regions where the 
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parallel flow has a large gradient the perp-parallel component R^n also enters. It is 
interesting to note that the direct dissipation in this model due to and Rs acting on 
the diamagnetic flow u* is exactly cancelled by part of the polarisation transfer, leaving 
only the time dependent and advective vorticity contributions, which constitute a real 
compressional effect. This does not reduce the computational complexity however. But 
since Eq. (73) says that V • Jp = 0, the polarisation current Jp may itself vanish in an 
enclosed domain, allowing the neglect of the work term involving Jp entirely. This only 
holds if the gradients of and Pi are set to zero at the outer boundary, however; if it 
is desired to incorporate a loss model at that boundary, then one must keep Jp • Vpi to 
conserve energy. The other unusual effects regarding the polarisation drift velocity are 
avoided because of the neglect of parallel transport in a one-dimensional model, but in a 
two-dimensional model they would have to be faced (see the comments preceding Eq. 36). 

It might be tempting to model strictly as flow generation, expecting it to balance 
the neoclassical friction in determining the equilibrium state of (f) [18]. But in many cases 
the only important effect on is R^ acting alone; there is evidence from experimental 
studies of internal transport barriers that in some regimes, especially coinjection of neutral 
beams, (the rotation or V-cross-B terms in the force balance are stronger than the pressure 
gradient and therefore) the Reynolds stress balances itself to zero in response to strong 
flows [19], and computations of drift wave turbulence in the presence of a background ExB 
vorticity show that the longer wavelengths are actually driven while the shorter ones are 
suppressed, so that very strong flow shear would provide a net drive for the turbulence and 
therefore self-deplete [20] . Such a self-regulator mechanism is probably best for modelling 
R^. Parallel flow shear is much less well studied, so the question of what to do for R_e|| 
is more open. 

Current transport codes neglect the partial time derivative in Eq. (73) and with it the 
difficulties in treating the polarisation transfer; the frictional damping then (or R^; itself) 
takes over the determination of ^ and the energy in uj^ is neglected. But in general can 
be as large as tty, so it is better to keep both flow energies. Actually, most transport codes 
neglect the question of plasma rotation entirely, as well as the evolution of the magnetic 
flux, which reduces the set of equations to Eqs. (70-72) for rig, Pe, and pi, solely with the 
particle and energy fluxes, which are easy to model consistently. Even in that case, we do 
still have the anomalous transfer terms due to the disturbances, and these could become 
interesting in several regions of the plasma. 

The existence of anomalous heat transfer has been noted before, in the context of 
a quasilinear model for the disturbances [21]. The more general version presented here 
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in the context of mean field theory incorporates models based upon linear instabilities as 
a subset, but it also allows incorporation of the results from general turbulence studies, 
including the strong mode structure changes brought about by the nonlinear dynamics 
of the ExB vorticity [7]. The quasilinear model indicates that Tgj is only important 
where T is significant; in turbulence however, especially electromagnetic turbulence, the 
robust activity in the dynamics involving Vy Jy ensures a prominent role for all three of 
the contributors to Tei in Eq. (79). Ultimately, further investigation and diagnosis of 
generalised turbulence computations will be required to validate any strong conclusions 
involving these efi^ects. 

V. The Characteristics of Transport Caused by ExB Turbulence 

In this section we confirm the main points of Section II concerning the special nature of 
ExB turbulent transport, compared with the more familiar random scattering via Coulomb 
collisions. 

The first point is the simplest: the ExB velocity is nearly divergence free, 

V • V£; !^ (82) 

and in a homogeneous magnetic field it is divergence free. For some laboratory configu- 
rations, cylindrical plasma geometry at very low beta bounded by plates at the ends, for 
example, this can be taken to be exact. The result is something well known from neu- 
tral fiuid dynamics: in a divergence free velocity field the pressure does not do work on 
the fiuid elements. The transfer of energy between the pressure and the velocity becomes 
small. For interchange turbulence the work done by the pressure on the ExB velocity 
divergence is responsible for the principal forcing effect on the velocity. But for the profiles 
this is overshadowed by the other transfer effects. A divergence free transporting velocity 
is something very different from the model in coUisional kinetics, by which the velocity is 
given by diffusion, 

UeVD = -D{ne,Te)Vne (83) 

for example. This velocity is only divergence free for specific dependences of D on the 
parameters. For example, the scale of the velocity divergence for a constant diffusion is 
the same as the scale of the density profile. But is exactly or nearly divergence free 
according to the magnetic geometry, regardless of what form ne(r) or Tg{r) take. 

The second point is the Poynting cancellation, Eq. (42). The presence of pVE in the 
transport with the factor of 3/2 is true for total energy, but what we want to know is if 
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that is just a result of formal manipulation. We recall that an ideal pressure equation may 
be written in terms of a transport term [(5/2)pv] and a work term [v- Vp], or alternatively 
in terms of a transport term [(3/2)pv] and a divergence [pV • v]. If the divergence is 
small, as we showed for v^;, then the second form is more germane. But it is even more 
interesting to know how the transport works for each piece of the energy rather than the 
total. The quantity (0J_l + pve) appears in the transport equation for the drift energy, 
Eq. 45. In modelling transport terms we usually write down the lowest order forms, which 
still conserve energy properly. We know that the diamagnetic current is at the same order 
as if TT-gcEj^ is comparable to Vp, so (f)J± is generally the same size as pvE- In fact 
they are close enough for the remnant to be negligible. A particularly useful form of the 
Poynting cancellation turns out to be the one with the diamagnetic current. 



This is a magnetic divergence, and it is now down to the next smallest row of terms in 
Eq. (50) as its scale is that of the magnetic field (i? in toroidal geometry), and the scale 
of interest is Lp. It is therefore a fundamental property of ExB flow dynamics that its 
contribution to the total energy transport is 



given that there are no radial equilibrium flows. The physical point of this is that the energy 
due to the fluid perpendicular velocities presents a transfer term of the form pV • v^, rather 
than the more usual form • Vp, to the thermal energy reservoir, due to the Poynting 
cancellation, so that the actual transport is due to the remaining {3/2)pvE- This role of 
the Poynting cancellation is what makes the ExB velocity somewhat "special" compared 
to all the others. 

The third point concerns magnetic flux diffusion. It is an interesting measured prop- 
erty of magnetically confined plasma experiments that the electron thermal transport is 
anomalous by more than two orders of magnitude, but the magnetic fiux diffusion is 
neoclassical [22] (classical diffusion, modified by toroidal drifts of particles whose parallel 
motion is in the long mean free path regime). If we write the mean field Ohm's law keeping 
magnetic fiutter, however. 



V-(0J*+pvs) = V- 



r c 



BxV(#) 



(84) 



Qe + Qi = 2 (p^e) 



(85) 




(86) 
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we instantly see why this should be true for ExB turbulence, If the magnetic disturbances 
are small or, as is more often the case, uncorrelated, then the anomalous flux diffusion is 
small; indeed it is zero for purely electrostatic turbulence. It is well known that (qe\\^±^ 
contributes negligibly to transport for drift wave turbulence in either slab or toroidal 
geometry, due to self consistency effects between and Tg [6,23]. Only when the current 
gradient is available as an energy source, as for tearing instabilities, [24], do these processes 
lead to appreciable transport. It then follows that the anomalous flux diffusion should be 
small. With determined by Ay, this holds as well for the current diffusion. Anomalous 
current diffusion can only result from the ExB nonlinearity in the electron inertia, and this 
is small compared to the mechanisms by which currents in a specifically current carrying 
plasma are generated. This is another nice contrast to classical diffusion, for in that case 
the flux diffusivity is much ^^reater than the particle diffusivity. Neglecting bj^, the mean 
fleld Ohm's law is purely diffusive, as we can see from 



,2 



£. + ^Vi^„ (87) 

having inserted Eq. (76) to make the flux diffusion explicit. Comparing the classical 
diffusivities (assuming a classical ryy), we flnd 

?7iic2 0.51 2 0.51 ^ 

Pe^e = (88) 



47r pe Pe 

With ~ 10~^ in the edge regions of fusion plasmas, this classical diffusivity is as large as 
the anomalous electron thermal diffusivity if not larger. An interesting consequence is that 
the MHD equilibrium in this regime is resistive. If the pressure gradient changes during a 
transport event it is not the case that the magnetic structure evolves isentropically through 
its successive equilibria. This point has been largely ignored in edge transport modelling, 
which takes the magnetic structure to be fixed. 

The fourth point is the various mechanisms of ExB fiow generation. Self consistent 
fiow generation by turbulence goes through the Reynolds stress [18], a fiuid dynamical 
mechanism. But for closed, toroidal geometry a sheared ExB velocity profile also means a 
radial electric field with radial dependence. In other words, a finite divergence and hence 
a finite charge density. But within low frequency fluid drift motion at scales larger than 
the Debye length, the charge density is always small, speciflcally, much less than UeC, the 
charge density of the electrons. Not only that, it is small in all phases of the dynamics. 
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which we can quickly show. Writing simple continuity equations for both fluids using the 
velocities in Eqs. (19,20), we find a charge continuity equation given by 

^ + V • (pchVE) + V • (J||b + J, + Jp) = (89) 

where pch = ngC — UiZe is the charge density. Taking into account that • V ~ d/dt, 
the first two terms are the same size. But the polarisation current can be written next to 
them, 

- V • {pch^rE) - V • Jp = V • (J||b + J,) (90) 

The polarisation divergence and the charge density term have the same form, and the 
polarisation divergence is larger by a factor of c^/f^, since Vj^0 = ^Trpch- All of the terms 
involving pch are then negligible, and we are left with V • J = as before. This underscores 
the strictness of the quasineutral character of the dynamics. It holds pointwise, not just 
in an averaged sense. We then conclude that any mechanism which should build up an 
electric field does so under the ExB fiuid dynamics and not due to transport of charges 
by currents. The total charge transport itself must vanish. This leads to what happens 
as a result of any charged particle source or sink mechanisms. It is quite possible to have 
source terms on the right side of Eq. (90), but instead of accumulating charge we then 
have a balance between sources and transport, 

-V-Jp = V- (J||b + J*) + >Se--5^ (91) 

so that the charge content remains small. This is the high throughput regime for electric 
charge, similar to the situation with free energy in drift wave turbulence where the current 
transfers thermal and drift free energy back and forth but the magnetic energy always 
remains small. In the high throughput regime for charges, quasineutrality is maintained 
and these source mechanisms act as a torque on the ExB vorticity [12]. Taking charges 
out of a flux tube causes the flux tube to rotate, since while the small but flnite V^(^ 
leads to no appreciable charge density, it the ExB vorticity it does lead to has a signiflcant 
role. Ultimately, the source of this angular momentum, as in any case for perpendicular 
momentum in fluid drift motion, is the background magnetic fleld. 

The last point concerns anomalous momentum transport. With the ExB vorticity 
proportional to the small but nonzero charge density, transport of angular momentum takes 
place through the Reynolds stress, which is nothing more than the nonlinear polarisation 
current. But in modelling the transport of this current one should not use an anomalous 
conductivity in a simple radial Ohm's law. Conductivity, or better its inverse, resistivity. 
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is a friction between the two fluids which conserves total momentum. Resistivity leads to 
particle diffusion, as the magnetic Lorentz force balancing resistive friction, 

vxB , , 

= r/i J (92) 

with perpendicular resistivity ryx = nieVel'^e^ ■, gives rise to a coUisional drift velocity 
which acts as a diffusion, once the MHD equilibrium constraint, JxB = cVp, is set in, 

VD = -^ViP (93) 

e 

(in the neglect of thermal forces). But in the momentum equation, the resistive friction 
cancels, as 

^nemev + V • (• • •) = Rje - ^Ve - n^e (^E + ^ xb) (94) 



dt 

add to become approximately 



— njMjU + V •(•••) = -Rie - Vpi + nj Ze (^E + ^ X b) (95) 



^pu + V •(•••) = -Vp + ^ (96) 

The only phenomenon left to transport momentum is the Reynolds stress (cf. Eq. 58), 
which is the dominant effect in the momentum fliix terms we did not explicitly write in 
Eqs. (94,95). The JxB force remains as a transfer mechanism between fluid and Poynting 
momentum, with the background magnetic held acting as an anchor. 

To summarise, transport of the background proflle quantities by small scale ExB 
turbulence can still be diffusive, provided the scale of motion is small compared to Lp. 
But its fundamental properties are qualitatively different from those of a kinetic diffusion 
via random thermal motions, and these should be taken into account when constructing 
transport models. An example model for one-dimensional transport including flows and 
magnetic induction has been given in Section IV. Attempts to extend conventional two- 
and three-dimensional edge transport models, which are able to treat the parallel fluxes as 
well as the radial transport since the proflles are no longer flux functions in their regime, 
to incorporate these effects are ongoing [25] . 
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